Kondo Lattice Scenario in Disordered Semiconductor Heterostructures 
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We study nuclear relaxation in the presence of localized electrons in a two-dimensional electron 
gas in a disordered delta-doped semiconductor heterostructure and show that this method can 
reliably probe their magnetic interactions and possible long-range order. In contrast, we argue that 
transport measurements, the commonly-employed tool, may not sometimes distinguish between 
spatial disorder and long-range order. We illustrate the utility of using the nuclear relaxation 
method to detect long-range order by analysing a recent proposal made on the basis of transport 
measurements, on the spontaneous formation of a two-dimensional Kondo lattice in a 2D electron 
gas in a heterostructure. 



Introduction- The possibility of long-range charge or 
magnetic order of strongly-correlated electrons in meso- 
scopic devices, such as Wigner crystals [3, [2], charge den- 
sity waves 0|, and Kondo lattices has attracted a 
great deal of attention in recent times. Theoretical inter- 
est in these systems stems from the low-dimensionality 
which enhances quantum effects, while the practical mo- 
tivation comes from the tunability of the material pa- 
rameters by electrical means, which is not achievable in 
bulk materials. Experimental probes for long-range or- 
der in mesoscopic devices have been usually based on 
transport measurements Q as their small size makes it 
difficult to employ standard bulk methods such as diffrac- 
tion and nuclear magnetic resonance (NMR){H, 0, 0|. 
However suitably adapted NMR methods are now be- 
ginning to emerge as very promising tools for study- 
ing electron interactions in mesoscopic systems - recent 
work shows that nuclear polarization may be generated 
[1, H, EH EH, E3| locally in such devices and its relaxation 
can be feasibly detected [H, ld| through two-terminal 
conductance measurements, and the behavior of the nu- 
clear relaxation rate conveys useful information about the 
electronic state in the device. For example in the context 
of the decade-old puzzle of the 0.7 conductance anomaly 
in quantum point-contact devices [H], NMR can be used 
to distinguish between three incompatible contesting the- 
ories - a Kondo effect, a spin-incoherent Luttinger liq- 
uid state, or a polarized electron liquid [l3| even though 
transport properties are similar in the three scenarios. 

In this paper, we study nuclear relaxation as a probe 
for long-range magnetic (and crystalline) order of local- 
ized spins in a disordered, metallic two-dimensional elec- 
tron gas (2DEG) in a delta-doped heterostructure. We 
find that the temperature dependence of the relaxation 
rate for a disordered few-impurity system approximately 
follows a linear— T law, while for strong enough inter-spin 
interactions, nuclear relaxation in a regular array, or a 
Kondo lattice, will show an exponential increase, e A / T , 
with decreasing temperature. In contrast, we argue that 
transport measurements will show no significant differ- 
ence between the two situations. As an application of our 
analysis, we discuss a recent experimental claim 0] based 



on transport measurements in disordered GaAs/AlGaAs 
delta-doped heterostructures, on the spontaneous forma- 
tion of a Kondo lattice in the 2D electron gas in the 
heterostructure. 

Experimental context- Kondo lattice materials, such as 
heavy fermion metals, are being intensely studied 3, 13] 
to understand the nature of the competition of the mag- 
netic ordering tendency of the localized electrons and the 
screening tendency (Kondo effect) of the conduction elec- 
trons, close to quantum criticality. A 2D Kondo lattice, 
if engineered in a heterotructure, would offer the twin 
advantages of reduced dimensionality and tunability of 
parameters [J], and, as we show below, nuclear relaxation 
can be used to study these systems. 

In Ref. Q it was observed that the 2DEG conductance 
showed an alternating splitting and merging of a zero 
bias anomaly (ZBA) upon varying the gate voltage V g . 
The authors interpreted these observations as evidence 
for the formation of a spin-1/2 Kondo lattice embedded 
in a 2DEG with the following physical picture. Vary- 
ing the gate voltage affects the 2DEG density, which, 
in turn, controls the sign of the RKKY exchange interac- 
tion, JnKKY{Rij) ~ (J 2 p/Rfj)cos(2kFRij), of the local- 
ized spins. Here J is the Kondo coupling of the localized 
spins with the conduction electrons, and p is the density 
of states at the Fermi energy. 

Nevertheless, the observation of the ZBA splitting is 
not sufficient to prove the existence of a Kondo lattice. 
Such an effect has been observed in the context of dou- 
ble quantum dot (DQD) systems 1|| and attributed 
to the competition of Kondo and interdot exchange in- 
teractions. Even in a sample with a small number of 
localized spins, the Kondo and RKKY competition will 
be dominated by the pairs of spins with the strongest 
exchange interactions. To this end, we need to show that 
the nuclear relaxation rates have qualitatively different 
signatures for the Kondo lattice and few Kondo impuri- 
ties scenarios. 

Nuclear relaxation takes place through nuclear cou- 
pling to localized spins S as well as conduction electrons 
o ■ Ti-xoc = A-dl ■ S + A S I ■ a. The relaxation contribu- 
tion from localized (electron) spins is usually much larger 
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Figure 1: (Color online) Plots showing the qualitative differ- 
ences in the temperature and (AFM) interdot exchange inter- 
action J ex dependencies of the nuclear relaxation rates T{~ for 
a double quantum dot system and Kondo lattice. Main plot: 
T\ (T) for (a) a double dot system; (b) a Kondo interaction 
dominated lattice (J ex /uiK < 1); (c,d) a Kondo lattice where 
J ex /ujK > 1 and T < (>)Tg { , where Tg lf is the mean-field 
transition temperature. Dotted curve interpolates between 
these two temperature regimes (there is no phase transition). 
Inset: T 1 _1 as a function of J ex /coK for (i) the double dot sys- 
tem - note that T^ 1 vanishes for J ex /ujK > tt; and (ii) for the 
Kondo lattice. 





FM 


AFM 


Double 
impurity 


Linear-T at low temp, 
and 1/T at high temp. 


Zero at low temp, and 
1/T at high temp. 


Lattice 


T/(T - T c ) 3/2 at high 
temp, and exp(l/T*) at 
low temp. 


T/{T - T c ) at high 
temp, and exp(l/T) 
at low temp. 



Table I: Distinguishing different physical scenarios with NMR 



in devices similar to those considered here [2l|. Taking 
only the localized spin part, the nuclear relaxation rate 
can be expressed in terms of the transverse impurity sus- 

ceptibihty, Tf 1 = J ^^\ m (^-)^. We show 

below that the temperature dependencies of 1/Ti in the 
Kondo lattice and few quantum dot scenarios are quali- 
tatively different. The results obtained in the paper are 
illustrated in Fig. Q] and Table U 

Model- We consider the following model Hamiltonian 
for S = 1/2 magnetic impurities Si in a 2-dimensional 
electron gas: 

H = J2^ c l« c ^ + J^i-Si- (1) 

k i 

Here <Ji the conduction electron spin density at r^. 

We use the "drone-fermion" representation for the lo- 
calized spins 0, IH: S+ = fJ X i/V2, S7 = X ifi/V2, 
and Sf = f}fi — 1/2, where fi and // are fermi- 
nonic operators and \ are rea l Majorana fermions de- 
fined by {xiiXj} — %• Note that the commutation 
relations for the impurity spins are automatically sat- 



isfied, obviating the need to impose local constraints 
on the fermion number. Introducing the bosonic op- 
erators, a l = (//cj T + XiCii/V2)/V2, h = (//c^ - 

Xicjj A/2) / V2, the interaction part of the Hamiltonian 
(up to constants and irrelevant terms) can be written as 

Hi„t = ~JJ2i ( a ! a i + ■ Factorizing H int using the 
Hubbard-Stratonovich transformation introducing fields 
A \ and A 2 , and further making the transformations 

Al,2 = |Al, 2 |e^.= , / - fe^~^\ c T (r) - c T (r)e*, 
c l( r ) — * c l( r ) e ^ 1 j the partition function can be written 
in path integral form: 

Z = j D(c, /, X , A) e" /o* dr[s 0+ s M+ ± E ( (|A*xl 2 +|A*| 2 )] . 
5 o = c i<r & + &) c k a + ^2(fid T fi + )pad T Xi), 

k<7 i 

S mt = l A ll i a i + 4) + l A 2l + + $ terms. (2) 

i 

We make a mean field analysis, neglecting the fluctu- 
ations in </>'s and A's. The frequency-dependent local 
transverse susceptib ility at low temperatures T <C Tk 
can be shown to be [23j 

xt KO = ( + (r) s - {Tl)) „ 2 — 

(ffsMBj TT (\uj m \ + LUK) 

Here oj m are bosonic Matsubara frequencies, and u>k — 
D exp(— 4/3/9 J). Note that ujk differs from the correct 
Kondo temperature, ksTx ~ De^ 1 ^ pJ \ This is an ar- 
tifact of the mean field approach. An analytic continu- 
ation, xf~(u m ) -> xt~{u) = {g s LiB) 2 h{-ifr^ + u) Ki ), 
to real frequencies leads to the well-known result T^ 1 — 
A 2 d k B T/TThL0 2 Kl for T<^lo k . 

Nuclei may also relax through their hyperfine coupling 
with conduction electrons. It is easy to see that at low 
temperatures, the ratio of the nuclear relaxation rates 
from impurity coupling and conduction electron coupling 
is (w^-p 2 i?g„7r) _1 , where R en is the electron-nucleus sep- 
aration. Thus the impurit y coupling mechanism domi- 
nates as long as R en < 1/ y ujkP^ 1 ^ 2 ~ 160nm. 

Double or a few-impurity system- We now consider 
the impurity spin susceptibility for two spins Si,S 2 
at R = Ri , R,2 which have an exchange interaction 
H ex = J ex (Ri2)Si ■ S 2 among them. If the wavefunc- 
tions of the localized electrons have a significant overlap, 
then direct exchange would be dominant. Indirect (or 
RKKY) exchange is more important at larger separa- 
tions. Here it also becomes important to compare the 
relative strengths of the RKKY interaction between the 
impurity spins with the hyperfine interaction of either 
of the impurities with neighboring nuclei. The RKKY 
interaction Jrkky falls off with distance R12 not faster 
than Jrkky ~ J 2 p/Ri2- This should be compared with 
Ad = A s /^ oc i?j ot , where i?dot is the size of the quan- 
tum dot in the plane of the heterostructure and Zi oc is 
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the thickness of the 2DEG. We use the following pa- 
rameters for a GaAs/AlGaAs heterostructure, Jp ~ 1, 
l ]oc ~ lnm, R dot ~ lOnm, A s = 3.8 x lCT 54 Jm 3 , and 
m = 0.063m e . Then Jrkky ^ A d /N nuc is satisfied if 
i?i2 <C y^l/AdP ~ 1mm, and this is true for most de- 
vices. Thus the nuclei couple to the RKKY bound pair 
Si + S2 rather than the spins separately. The impurity 
susceptibility now involves both on-site and intersite cor- 
relations, and we have, to leading order in inter-impurity 
interaction for T <C Tk, 



rr A = 



A 2 d k B T 
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(4) 



When J ex /tOK > tt, the nuclear relaxation rate is sup- 
pressed to zero: this is the maximum value of J ex /u>K for 
which the behavior is governed by the Kondo screening 
of the impurity spins. Indeed, even for a large ferromag- 
netic coupling of the spins, the ground state is a Kondo 
singlet [2J|. At antiferromagnetic couplings J ex > nujK, 
the ground state is an RKKY singlet which is unable to 
exchange spins with the nuclei. While our analysis is only 



to leading order in J ex , more accurate calculations .24] 
based on numerical renormalization group methods have 
shown that this critical point occurs at J ex /k B TK ~ 2.2. 

We discuss now the validity of the mean field treat- 
ment. Note that the mean field A corresponds to bind- 
ing energy of the impurity fermion with the local con- 
duction electron. Ignoring fluctuations of the phase of 
A results in underestimation of ujk (see text following 
Eq. [3] and also Ref. [13]). Amplitude fluctuations of 
A may be ignored as long as the Kondo energy domi- 
nates inter-impurity exchange, i.e., for small J ex /u>K- In 
fact, the mean field approach is incapable of capturing 
the physics of the magnetically ordered phase. 

With a larger number of spatially disordered impurity 
spins, one can show that for weak inter-impurity inter- 
actions, the nuclear relaxation rates have linear-T be- 
havior with logarithmic factors arising from the random 
distribution of Kondo temperatures of individual impu- 
rities 25|. For strong inter- impurity exchange interac- 
tions, we can ignore the Kondo effect to leading order. In 
that case, it is known that the magnetic susceptibility at 
low temperatures is dominated by pairs with the weak- 
est exchange interactions [20j - this leads to a weakly- 
increasing susceptibility e cin 1 ( T °/ T ) (instead of zero for 
the double impurity case). Nevertheless, the nuclear re- 
laxation rate is dominated by the linear-T prefactor as 
the exponential term is weaker than any power law. 

Kondo lattice- The main physical difference from the 
two-impurity case is the existence of low energy mag- 
netic excitations in the lattice for any value of the ratio 
J ex /ujK- As a result, significant nuclear relaxation still oc- 
curs for large antiferromagnetic inter-impurity couplings 
unlike the two-impurity case where it vanishes. 

We consider first the scenario where we have a lattice 
of Kondo impurities with a weak exchange interaction 



( J ex <C ujk ) among the neighboring spins. Suppose that 
J ex (q) has maximum value at q = Q, and assume the 
wavevector dependence in the vicinity of the maximum 
is J ex (Q + q) = J ex (Q) - (D s /n- lmp )a 2 q 2 , where a is the 
lattice constant of the Kondo array and D s the spin wave 
stiffness. The random phase approximation (RPA) sus- 
ceptibility in this momentum region has the form 



(.dstlBf 



7T (uj sf (T) 



ifku 



D s a 2 q 2 



(5) 



where u>. f (T) = 



Jex(Q)"in 



A new energy 



scale w s /(0) = ujk — Jex(Q) n imp I ' 7r appears representing 
the competition of Kondo and inter-impurity exchange 



interactions. As J ex (Q)f 



nuiK, the uniform, static 



transverse susceptibility tends to diverge signaling a mag- 
netic phase transition. Using the known temperature 
dependence of the susceptibility of a Kondo impurity, 
xf~(0,T) ~ Xi(0)(l - Ck 2 B T 2 /uj 2 K ),(C is a constant of 
order 1), together with the frequency dependence of Im% 
from Eq. \E[ the nuclear relaxation rate for k B T <C luk 
turns out to be 



Tf 1 = A 2 d k B T/Aii 2 hD 8 Lu s f{T). 



(6) 



There is a crucial difference between the nuclear relax- 
ation results for the Kondo lattice in Eq. [6] and the two- 
impurity case. Consider for simplicity ujki — lok2 = ^k- 
First, near the transition J ex {Q)n\ mv /ujK = 7r, 1/Ti for 
the Kondo lattice is large and finite, while it tends to 
vanish for the two-impurity case. 

Now we consider the case when localized spin-spin 
interaction is dominant and neglect the Kondo inter- 
action in the zeroth order. We are particularly inter- 
ested in the regime close to a magnetic phase transition. 
The Hamiltonian describing the system would be H = 
Efe, CT Cfc4<T C ^ + E q J ex(q)S q -S_ q + JX; i o- l -S i , where J 
is to be treated now as a perturbation. J ex (q) represents 
all exchange processes except indirect exchange (RKKY), 

jRK KY (q,-) = J 2 E k ^::::::;: +lS - ^ we m ay 

write the effective inter-impurity exchange interaction as 
J ex (<2,cj) = J ex (<7) + Jrkky(<7,^)- where We now Tay- 
lor expand the exchange interaction near its extremum, 
Jex(Q + q) - Jex(Q)(l ~ a 2 q 2 ), where J ex (Q)a 2 = 
(D s /rii mp )a 2 . The spin susceptibility near the ordering 
point is approximately 



Ak B T™ { {a 2 /e + a 2 q 2 - z7Q+ q M} 



-(7) 



where Q is the wave vector of ordering, T™ f is the mean- 
field magnetic transition temperature, £(T) is the mag- 
netic correlation length and 7 q (w) is the imaginary part 
of Jex(q,w)/J ex (Q), 7qM - ^(Jp) 2 fkj/4J ex (Q)k F q = 
7(q)w. We can now estimate the nuclear relaxation rate. 
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For an antiferromagnetic (AFM) square lattice, the or- 
dering happens at Q = (x/a, 7r/a). Then the relaxation 
rate at the site of any given impurity is 



T 1 -\T) = 



A 2 d (J P ) 2 JUQhf mp T £(T) 



UhD 2 k F Q 



T"mf 



(8) 



Eq. [8] differs from estimates [2g] of 1 /Ti for Heisenberg 
antiferromagnets because in our case, the magnon de- 
cay is on account of the RKKY coupling of the impurity 
spins. Similarly for the ferromagnetic (FM) case, 



T -i (T) 



A 2 d n(Jp) 2 Jex(<9)n irap T 



(9) 



l2mk F aD 2 s 

The temperature dependencies of the correlation lengths 
are similar for the AFM and FM cases, 



\ aexp(27rD^/fc B T), 



^T 1 -v. a^mf 

J > J c . ( 10 ) 



where the low temperature behavior for the antiferro- 
magnet was obtained in Refs. 2fL 13], and for the ferro- 
magnet from Refs. [13, 01- D' s « 0.18J ex (Q)n imp is the 



exact spin wave stiffness at T = for a 2D (square lat- 
tice) Heisenberg magnet. These results also differ from 
the usually-encountered 3D Kondo lattice systems [29| . 
because of the qualitative difference in the behavior of 
£(T) at low temperatures. 

Finally let us discuss the effect of the Kondo inter- 
action on our results. In presence of inter-impurity 
exchange interactions, the singular Kondo corrections 
(~ {J p)\n{D / ksT)) to the gyromagnetic ratio of the 
impurity spins are modified to (Jp) ln(.D/yjJ + k 2 B T 2 ) 
[30|| . Consequently, the primary effect of Kondo correc- 
tions is to decrease the Stoner critical temperature T™ f 
as well as the pre-factor in the expressions for the nuclear 
relaxation rates but the temperature dependence of T^ 1 
does not change significantly. Eqs. [4j El [8] and [9] are our 
main results and are plotted in Fig. [TJ 

There is also a possibility of a magnetic instability be- 
low the Kondo temperature. In this case the ZBA split- 
ting would be smaller than the "heavy fermion" band- 
width, wjc, unlike the above case where the ZBA splitting 
is larger. As an example, for a Fermi liquid with FM spin 
fluctuations, results for nuclear relaxation available in the 



literature [3lJ| and are similar to our Jrkky / <^K > 1 re- 
sults. The difference will be seen in the magnitude of 
ZBA splitting relative to luk ■ 

In summary, we calculated the nuclear relaxation rates 
T^ 1 for the Kondo lattice and the few disordered mag- 
netic impurities cases and showed that they have qualita- 
tively different low temperature behaviors: when inter- 
spin exchange interactions are strong compared to the 
Kondo energy luk, the temperature dependence of T^ 1 
for the few-impurity system will follow an approximate 
linear— T law, while for the Kondo lattice T^ 1 will show 



an exponential behviour e A l T at low temperatures. In 
contrast, we argued that transport measurements 0| in 
this case may not provide a clinching evidence for the 
formation of crystalline order (Kondo lattice) . The expo- 
nential temperature dependence is special to two dimen- 
sions and indicates stronger spin fluctuations: a power- 
law behavior is expected in three dimensions on either 
side of the transition temperature[29]. These results also 
differ from a 2D Heisenberg magnet because in our case, 
magnon decay is mediated by conduction electrons. We 
hope our study will work towards encouraging the use of 
NMR measurements as an additional handle for study- 
ing magnetism and long-range order in low-dimensional 
conductors. 
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